Electromechanical Effects in Carbon Nanotubes 
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We perform ab initio calculations of charged graphene and single- wall carbon nanotubes (CNTs). 
A wealth of electromechanical behaviors is obtained: (1) Both nanotubes and graphene expand 
upon electron injection. (2) Upon hole injection, metallic nanotubes and graphene display a non- 
monotonic behavior: Upon increasing hole densities, the lattice constant initially contracts, reaches 
a minimum, and then starts to expand. The hole densities at minimum lattice constants are 0.3 
e|/atom for graphene and between 0.1 and 0.3 |e|/atom for the metallic nanotubes studied. (3) 
Semiconducting CNTs with small diameters (d < 20 A) always expand upon hole injection; (4) 
Semiconducting CNTs with large diameters (d > 20 A) display a behavior intermediate between 
those of metallic and large-gap CNTs. (5) The strain versus extra charge displays a linear plus power- 
law behavior, with characteristic exponents for graphene, metallic, and semiconducting CNTs. All 
these features are physically understood within a simple tight-binding total-energy model. 

PACS numbers: PACS numbers: 71.20.-b, 73.22.-f, 73.90.+f 



More than 10 years after their discovery by Iijima JjJ, 
carbon nanotubes (CNTs) are still capable of unveiling 
surprising and fascinating physical properties. An impor- 
tant feature of CNTs is the dependence of their electronic 
properties on diameter and chirality ||, ||, ||, |], roj . How- 
ever, not only the electronic properties of CNTs are at- 
tractive, as these materials have superb mechanical prop- 
erties as well |t| , making their use in nanomechanical ap- 
plications very promising Q . The possibility of electronic 
control of nanomechanical devices would certainly bring 
a great improvement in their effective switching times, 
control and precision H. 

For all these reasons, it is important to study the elec- 
tromechanical properties of CNTs. A few experimental 
works have focused on electrostatically-driven mechan- 
ical responses of CNTs 0, 0. However, mechanical 
response of CNTs can be driven not only by electrostat- 
ics, but also by quantum-mechanical effects. Indeed, it 
is well known that the in-plane lattice constant of in- 
tercalated graphite expands or contracts relative to pure 
graphite because of charge-transfer effects [|l2[ Re- 
cently, the mechanical response of "nanotube sheets" - 
sheets of entangled single- wall CNTs bundles - to electro- 
chemical charge injection has been investigated, and its 
use in actuators has been proposed E3] . In that work, the 
electromechanical response of CNTs has been found to be 
stronger and more non-linear than that of graphite. Con- 
trolling and optimizing the electromechanical response 
of nanotube sheets is a difficult task, since it involves 
problems like bundling and entanglement of CNTs. For 
these reasons, it would be highly desirable to investigate 
the electromechanical properties of individual nanotubes 
and their dependence on chirality and diameter. Since 
experimental control of these parameters is still beyond 
the current state-of-art, the predictive power of ab initio 



theoretical calculations makes them the method of choice 
for such studies. 

Initially, we perform ab initio calculations of charged 
graphene and single-wall metallic and semiconducting 
CNTs. Our calculations are performed within the 
Density Functional Theory (DFT) and pseudopotential 
frameworks with a numerical-atomic-orbitals basis set 
(SIESTA code) Jl|. SIESTA has been successfully used 
in a number of studies on CNTs jl6| [l7) . 

Fig. |l|(a) shows the relative variation of the in- 
plane lattice parameter for graphene or the unit cell axial 
length for metallic zig-zag (12,0) and armchair (5,5) 
CNTs, as a function of the extra charge |L9| (in units 
of e|) per atom, q. Negative q (extra electrons) cause 
expansion of graphene and metallic CNTs, whereas low 
positive q (holes) cause contraction. Large amounts of 
hole injection eventually cause expansion of graphene and 
CNTs, in qualitative agreement with experiments [ pd| . 
Fig. |l](a) also shows, from the comparison between the 
results for the (5,5) and the (12,0) CNTs, that the elec- 
tromechanical response of metallic CNTs is sensitive to 
chirality and/or diameter. Also in agreement with exper- 
imental observations, the strain ^ vs. q appears to be 
more non-linear for CNTs than for graphene. The largest 
negative strain is -0.05% for graphene and -0.06% for 
(12,0) nanotubes, therefore larger than the experimen- 
tal value of -0.02 to -0.03%, indicating that the effect of 
bundling in nanotube sheets may hinder the electrome- 
chanical response of CNTs. A linear fit in the q < 
regime for graphene gives j£- = 0.060, in good agree- 
ment with the experimental value of 0.066 for interca- 
lated graphite 0. 

Fig. |l|(b) shows the behavior of j~ vs. q for a semi- 
conducting (11,0) CNT. Surprisingly, the CNT expands 
regardless of the sign of q, in contrast to what is observed 
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bands of undistorted graphene are 



£o(fc) = £2 P ±t f(k), 



(1) 
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FIG. 1: (a) Relative variation of the lattice constant (for 
graphene) or the axial unit cell lengths (for metallic CNTs), 
jji, as a function of the extra charge, q. Points are ab initio 
results, lines are best fits using Eqs. @ and (|l|). Jb) Same 
as (a) for semiconducting CNTs and fit using Eq. (D_4). Inset: 
plot of ^-{q = 0.01) for semiconducting CNTs of different 
diameters, showing the crossover to "metallic-like" behavior, 
(c) Log-log plots showing the different power laws of 8L(q), as 
discussed in the text. The slopes are 0.52 ± 0.04, 0.94 ± 0.15, 
1.32 ± 0.21 and 1.90 ± 0.15, for graphene, (12,0), (5,5) and 
(11,0) CNTs respectively, consistent with the predictions of 
the TB model ((a - 1), where a is defined in Eq. @): 1/2, 
1, 1 and 2. 



for graphene and metallic CNTs. 

A key element to understand these results is the en- 
ergy positioning of bonding and anti-bonding states in 
graphene and CNTs. A widely accepted view, based on 
a 7r-orbital, nearest-neighbor (nn) tight-binding (TB) de- 
scription of the energy bands of graphene Q , is that the 
lower (valence or ir) band has a bonding character and 
the upper (conduction or 7r*) band has an anti-bonding 
character. Indeed, within such a description, the energy 



where e 2p is an "atomic level" for the 2p state of the C 
atom (corrected by the crystalline potential), to is the 
nn hopping integral for an undistorted graphene lattice 
(lattice constant ao), and f(k) is a function of the wave 
vector in the 2-d Brillouin Zone (BZ) |Q. In this case, e 2p 
coincides with the Fermi level, Ep, for a neutral system 
at the K points of the BZ, where f(k) = 0. Therefore, 
extra electrons would occupy anti-bonding states and ex- 
tra holes would go to bonding states. In cither case the 
lattice would expand. Therefore, the nn TB approxi- 
mation cannot explain the contraction of graphene and 
metallic CNTs for low positive q. 

However, as pointed out by Kertesz |2^] , inclusion of 
next-nearest-neighbor (nnn) interactions lifts the energy 
eigenvalues at the K points above e 2p , therefore making 
these states antibonding. Since our method is based on 
an atomic-orbital basis, we can provide ab initio verifi- 
cation of these ideas. From the ab initio Hamiltonian 
(with interactions up to 4 th nn), we estimate e 2p to be 
roughly 1.6 eV below Ep for graphene. This is illustrated 
in Fig. ||. Therefore, the states around Ep are antibond- 
ing for both bands, and now if electrons are gradually 
removed from the n band (positive charging), the bond 
lengths will contract, in agreement with experiments |l~^ ] 
and with the results in Fig. [I]. If positive charging con- 
tinues, electrons will be eventually removed from states 
below e 2p (therefore from bonding states) and the bond 
lengths will start to increase again, also in agreement 
with the ab initio results. 

We explain this behavior in a more quantitative way 
via a total-energy TB model. The total energy of 
graphene or CNTs is a function of the number of elec- 
trons per atom, N, and the bond length variation 5L, 
and it is expressed as 



E? ot (5L) =E b + E r =2j2 e^(fc) + E r 



(2) 



where Eb is the band energy and E r is a repulsive energy. 
The summation runs over occupied states (accounting for 
the proper value of N) and ssL^k) is the band structure 
for a deformation SL . 

From Eq. (Q), we can write the band structure depen- 
dence on the deformation as 



£SL{k) = e 2p + {e (k) - e 2p ) 



To' 1 



(3) 



where t$p is the hopping matrix element for a distorted 
lattice. Assuming a power-law dependence of t on the 
bond length L, i.e. (t$L — to(Lo/ L) a ) and taking the 
limit of small 8L, we can linearize: 



£sh{k) = e (k) - (5L(e (k) - e 2p ), 



(4) 
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where ( — a/Lo. Let us consider the case of a neutral 
system, with iVo electrons. From Eqs. @ and (||), Ef, is 
given by 



sJAkf 



E^(SL) = E»°(0) C6L[E?°(0) N e 2p }. (5) 

For a neutral system, by definition, the total energy must 
be minimum at SL = 0, so 



E*%(8L) = E%(0)+j3{6L)* 



(6) 



where (i is an elastic constant. Substitution of (|5|) and 
@ in (||) yields an expression for the repulsive energy of 
a neutral system at arbitrary SL: 

E?°(6L) = E^(0) + P(5L) 2 + (5L(E^(0)-N e 2p ) (7) 

Consider now a system with an arbitrary number of 
electrons N. The calculation of E^ is similar to that for 
the neutral case, yielding 



E» (SL) = E» (0) - C6L(E? (0) - Ne 2p ) 



(8) 



We also assume that, to first order, the dependence of 
E r on the bond length is not very much affected by the 
extra charge. Then, from Eq. (Q) we have 

E?(5L) = E?(0) + (3(SL) 2 + CSL(E b N °(0) - N e 2p ) (9) 

Combining Eqs. (|), (|) and (§, and defining the extra 
charge as q = —AN = No — N, wc finally arrive at an 
expression for the total energy of a charged system: 

El t (5L) = ElM + 0{6Lf + (6Lq(AE b - e 2p ). (10) 



Here, AE h 



<(0)-<°(0) 
AJV 



is the variation in band en- 



ergy per extra electron. By imposing gj: 



0, 



obtain the unit cell length distortion as a function of q: 



SL(q) 



2f3 l 



AE h 



£2p 



(11) 



Eq. ( pd| ) has a very simple physical meaning. Consider 
for instance electron injection (q < 0). If AEb — e 2p > 0, 
extra electrons are predominantly placed in anti-bonding 
levels, so the lattice expands (SL > 0). Similar argu- 
ments can be used for electrons in bonding levels (yield- 
ing lattice contraction) and holes in anti-bonding or 
bonding levels (yielding contraction or expansion). 

We now calculate SL(q) for graphene and CNTs. Let 
us consider graphene first. In the vicinity of Ep, we 
can use a simple "conical" approximation for the energy 
bands, e(k) — Ep ± j\k\, where -I- stands for electrons 
and — for holes, and the origin in fc-space is shifted to the 
K point. AEb may then be obtained by integrating e(k) 
over the region of the BZ occupied with extra electrons, 
i.e. a disk of radius Afc which is related to q as ^fq = 



, where s is the area per atom. After integrating 



and substituting the result into Eq. (|ll|), we get 



(12) 



The energy bands of CNTs are obtained by slicing the 
energy bands of graphene along the proper quantization 
lines 0. We show schematically in Fig. || the result of 
this procedure in the vicinity of one of the K points. For 
metallic CNTs, one of the quantization lines will cross 
the K point, yielding linear 1-d bands. If the CNT is 
semiconducting, the 1-d bands will be parabolic. 

Semiconductor 



Mclal 




FIG. 2: (color) Schematic band structure of graphene and 
CNTs, in the vicinity of a if point. The green and red por- 
tions of the energy surface indicate bonding and anti-bonding 
states (energy below and above £2p), respectively. The two 
vertical planes indicate the effects of fc-space quantization for 
metallic and semiconducting CNTs. 

For metallic CNTs, the calculation of AEf, is very sim- 
ilar, the only essential difference being that the integra- 
tion is now one-dimensional and the integration limits 
are related to the q as =fq = where £ is the CNT 
length per atom. The result is 



c 



rrrq\ 



8L(q) = --q[E F -e 2p ^ J 



(13) 



Finally, for semiconducting CNTs we use a parabolic 
dispersion e(k) = Ef ± ^ ± 77/c 2 , where e g is the energy 
gap. The result for 8L(q) is 



E F ±^-e 2p ± gL, 



(14) 



From Eqs. @, Q and @) we see that the elec- 
tromechanical behaviors of graphene and CNTs can all 
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be described by simple expressions: 



SL(q) 



= -aq + b\q\ c 



(15) 



where the linear coefficient a controls the electromechan- 
ical response for low charge injection and depends pre- 
cisely on the positioning of Ef with respect to €2 P - The 
non-linear term is defined by the exponent a, which is 
3/2, 2 and 3 for graphene, metallic and semiconducting 
CNTs. We test these predictions from the ab initio re- 
sults. In Fig. |l|(c) we display log-log plots of j^- + a 
vs. \q\ for each system. Within error bars, the exponents 
predicted by the simple TB model are beautifully con- 
firmed by the ab initio results. Notice that both (5,5) and 
(12,0) metallic CNTs have the same exponent a = 2. The 
stronger non-linearity in strain vs. q for the (5,5) CNT - 
evident from Fig. |l|(a) - is due to a higher coefficient of 
the quadratic term in j^{q)- Such a large difference in 
the quadratic term indicates a strong dependence of the 
non-linear response of metallic CNTs on their chirality. 

Let us consider again semiconducting CNTs. From 
Fig. H and Eq. ([Til), we see that, for large enough e g , the 
upper valence states will stay below ti v . Therefore, for 
positive q, electrons will be removed from bonding states 
and expansion will occur. This is precisely what Fig. |l](b) 
shows. On the other hand, we can predict that if e g is 
small enough, semiconducting CNTs will behave qual- 
itatively like metallic ones, i.e., contracting upon hole 
injection. Since, e g oc 1/d 2 (where d is the CNT diam- 
eter) JtJ, we test this prediction by performing ab initio 
calculations on zig-zag semiconducting CNTs of differ- 
ent diameters. The inset of Fig. [l](b) shows a plot of 
SL/Lq as a function of d for a fixed, small extra charge 
of 0.01 holes/atom. This indeed shows that the decrease 
in the gap of (n, 0) semiconducting CNTs with increas- 
ing n ultimately leads to a recovery of metallic nanotube 
behavior. 

In conclusion, we have performed ab initio studies of 
electromechanical effects in graphene and CNTs. The 
diversity on electromechanical responses is fascinating, 
and it has an amazingly simple explanation in terms 
of the bonding and anti-bonding nature of the elec- 
tronic states within a TB model. The observed stronger 
non-linear behavior in SL(q) for CNTs as compared to 
graphene is reproduced and explained, arising naturally 
as a dimensionality effect in the distribution of extra elec- 
trons over the BZ's of graphene and CNTs. Our results 
provide important guidelines for fabrication of nano- 
electromechanical devices. Indeed, when experimental 
control of production and assembling of nanotubes with 
specific chiralities is achieved, the possibilities in the ar- 
chitecture of different devices by combining the variety 
of effects described here will be only limited by our imag- 
ination. 
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Note added: After the submission of this work, a re- 



lated study by Gartstein et al 24 has appeared in the 
literature. 
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